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A Hamiltonian is said to be quasi- exactly solvable (QES) if some of the energy 
levels and the corresponding eigenfunctions can be calculated exactly and in closed 
form. An entirely new class of QES Hamiltonians having sextic polynomial potentials 
is constructed. These new Hamiltonians are different from the sextic QES Hamilto- 
nians in the literature because their eigenfunctions obey "PT-symmetric rather than 
. Hermitian boundary conditions. These new Hamiltonians present a novel problem 

£^ \ that is not encountered when the Hamiltonian is Hermitian: It is necessary to distin- 

guish between the parametric region of unbroken VT symmetry, in which all of the 
eigenvalues are real, and the region of broken VT symmetry, in which some of the 
eigenvalues are complex. The precise location of the boundary between these two 
regions is determined numerically using extrapolation techniques and analytically 
£> ' using WKB analysis. 
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IT) ■ I. SEXTIC QES HAMILTONIANS 

o . 

The purpose of this paper is to introduce a new class of quasi-exactly solvable (QES) 
Hamiltonians having sextic polynomial potentials. While these new kinds of QES Hamil- 
tonians have positive, real eigenvalues, they have not yet been discussed in the literature 
because they are not Hermitian. Instead, they are VT symmetric. 
q-i The term quasi-exactly solvable (QES) is used to describe a quantum-mechanical Hamil- 

tonian when a finite portion of its energy spectrum and associated eigenfunctions can be 
found exactly and in closed form pj. Typically, QES potentials depend on a parameter J, 
and for positive integer values of J one can find exactly the first J eigenvalues and eigen- 
functions, usually of a given parity. It has been shown that QES systems can be classified by 
using an algebraic approach in which the Hamiltonian is expressed in terms of the generators 
of a Lie algebra 0,131113. 

Perhaps the simplest example of a QES Hamiltonian having a sextic potential is 0, 0] 

H = p 2 + x 6 -(47 -l)x 2 , (1) 

where J is a positive integer. For each positive integer value of J, the time-independent 
Schrodinger equation for this Hamiltonian, 

-ip"(x) + [x 6 - (47 - l)x 2 ]^(x) = Eip(x), (2) 

has J even-parity eigenfunctions in the form of an exponential times a polynomial: 

ij(x)=e-^'j2c k x 2k . (3) 

k=0 
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The polynomial coefficients c k (0 < k < J — 1) satisfy the recursion relation 

4(J - k)c k ^ + Ec k + 2{k + l){2k + l)c fc+ i = 0, (4) 

where we define c_i = cj = 0. The simultaneous linear equations (J3J) have a nontrivial 
solution for Co, Ci, cj_i if the determinant of the coefficients vanishes. This determinant 
is a polynomial of degree J in the variable E. The roots of this polynomial are all real 
and are the J quasi-exact energy eigenvalues of the Hamiltonian Note that all of the 
QES eigenfunctions (p} ]) of H in ([T jl have the form of a decaying exponential exp(— \x A ) 
multiplying a polynomial. This is the standard form in the literature for the eigenfunctions 
of any QES Hamiltonian whose potential is a polynomial. 

The QES Hamiltonians associated with Hermitian Hamiltonians have been examined in 
depth and classified exhaustively pj. However, in 1998 new kinds of Hamiltonians that 
have positive real energy levels were discovered [H H|. These new kinds of Hamiltonians 
are not Hermitian (H ^ W) in the usual Dirac sense, where the Dirac adjoint symbol 
f represents combined transpose and complex conjugation. Instead, these Hamiltonians 
possess VT symmetry H = H vr ; that is, they remain invariant under combined space and 
time reflection. This new class of non-Hermitian Hamiltonians has been studied heavily 0, 
11, l3| and it has been shown that when the V T sy mmetry is not broken, such Hamiltonians 
define unitary theories of quantum mechanics 3 

The key difference between Hermitian Hamiltonians and complex, non-Hermitian, VT- 
symmetric Hamiltonians is that with PT-symmetric Hamiltonians the boundary conditions 
on the eigenfunctions (the solutions to the time-independent Schrodinger equation) are im- 
posed in wedges in the complex plane. Sometimes these wedges do not include the real 
axis. (A detailed discussion of the complex asymptotic behavior of solutions to eigenvalue 
problems may be found in Ref. 



The discovery of "PT-symmetric Hamiltonians was followed immediately by the discovery 
of a new class of QES models. Until 1998 it was thought that if the potential associated 
with a QES Hamiltonian was a polynomial, then this polynomial had to be at least sextic; 
its degree could not be less than six. This property is in fact true for Hamiltonians that 
are Hermitian. However, in 1998 it was discovered that it is possible to have a QES non- 
Hermitian complex Hamiltonian whose potential is quartic |l5| : 

H = p 2 — x 4 + 2iax 3 + (a 2 - 2b)x 2 + 2i(ab - J)x. (5) 

Here, a and b are real parameters and J is a positive integer. For a large region of the 
parameters a and b, the energy levels of this family of quartic Hamiltonians are real, discrete, 
and bounded below, and the quasi-exact portion of the spectra consists of the lowest J 
eigenvalues. Like the eigenvalues of the Hamiltonian the lowest J eigenvalues of these 
potentials are the roots of a Jth-degree polynomial [if) . 

The reality of the eigenvalues of H in (JHJ) is ensured by the boundary conditions that its 
eigenfunctions are required to satisfy. The eigenfunctions are required to vanish as \x\ — ■> oo 
in the complex-x plane inside of two wedges called Stokes wedges. The right wedge is 
bounded above and below by lines at 0° and —60° and the left wedge is bounded above and 
below by lines at —180° and —120°. The leading asymptotic behavior of the wave function 
inside these wedges is given by 

il){x) ~ e"" 3/3 (|ar| -> oo). (6) 
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The new class of QES sextic Hamiltonians reported in this paper has the form 

H = p 2 + X 6 + 2ax A + (47 - 1 + a 2 )x 2 , (7) 

where J is a positive integer and a is a real parameter. These Hamiltonians are very similar 
in structure to those in $IJ) and to the other QES sextic Hamiltonians discussed in the 
literature [l| , but their distinguishing characteristic is that the asymptotic behavior of their 
eigenfunctions in the complex-x plane is different. 

Let us examine first the asymptotic behavior of the eigenfunction solutions to the 
Schrodinger equation (J2J). For brevity, we call the eigenfunctions in (J3J) the good solutions to 
(J2J) because they satisfy the physical requirement of being quadratically integrable. These 
good solutions decay exponentially like exp(— \x A ) as x —>■ ±oo, while the corresponding 
linearly independent bad solutions grow exponentially like exp(|x 4 ) as x — > ±oo. In the 
complex-x plane the good solutions (JHJ) decay exponentially as \x\ — > oo in two Stokes 
wedges that are centered about the positive and the negative real-x axes. These wedges 
have an angular opening of 45°. The bad solutions grow exponentially in these wedges. At 
the upper and lower edges of these wedges the good and bad solutions cease to decay and 
to grow exponentially and they become purely oscillatory. 

As we move downward past the lower edges of these wedges, we enter a new pair of 
Stokes wedges. These wedges also have a 45° angular opening and are centered about the 
lines argx = —45° and argx = —135°. In these lower wedges, the good solutions grow 
exponentially as \x\ — > oo and thus they behave like a bad solutions. 

In the lower pair of wedges we can find solutions to the new class of Hamiltonians in (J7J) 
that behave like good solutions. These new PT-symmetric eigenfunctions have the general 
form of the exponential exp(|x 4 + \ax 2 ) multiplied by a polynomial |l7j ]: 

j-i 

<4,(x) = e xA ' A+ax2 l 2 Y,CkX 2k . (8) 

fc=0 

Hamiltonians having even sextic polynomial potentials are special because such Hamil- 
tonians can be either Hermitian or PT-symmetric depending on whether the eigenfunctions 
are required to vanish exponentially in the 45° wedges containing the positive and nega- 
tive real-x axes or in the other pair of 45° wedges contiguous to and lying just below these 
wedges in the complex-a; plane. The solutions for these two different boundary conditions 
are somewhat related. Specifically, a good solution in one pair of wedges becomes a bad 
solution in the other pair of wedges. However, a bad solution in one pair of wedges does not 
become a good solution in the other pair of wedges, as we now explain. 

Given a g ood solution ip goo &(x) in one pair of wedges, we use the method of reduction 
of order |l8( to find the bad solution. We seek a bad solution in the form iphad{ x ) = 
ipgood{ x ) u ( x ) > where u(x) is an unknown function to be determined. Substituting the bad 
solution into the Schrodinger equation —tp"(x) + V(x)ip(x) = Eijj(x), we get the differential 
equation satisfied by u(x): 

ip goo d(x)u"(x) + 2tp' good (x)u' \x) = 0. (9) 

We solve this equation by multiplying by the integrating factor Vgood(^) and obtain the 
result 

/ r x 

-2 



Vw(z) = VWi(z) / ds [^good(s)] +C), (10) 
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where C is an arbitrary constant. 

This bad solution always grows exponentially in the two wedges in which the good solution 
decays exponentially. How does this bad solution behave in the other pair of wedges in which 
the good solution grows exponentially? We can always choose the constant C so that the 
bad solution vanishes as \x\ — > oo in one of these two wedges. However, in the other of the 
two wedges, the bad solution will always grow exponentially. Thus, while the good solution 
becomes bad as we cross from one pair of wedges to the other, the bad solution does not 
become good. 

II. DETERMINATION OF THE VT BOUNDARY 

The difference between the Hermitian Hamiltonians in JTJ and the PT-symmetric Hamil- 
tonians in J7J) is that the Hermitian Hamiltonians always have real eigenvalues. The VT- 
symmetric Hamiltonians in (J7J) have real eigenvalues only if the VT symmetry is unbroken; 
if the VT symmetry is broken, some of the eigenvalues will be complex. Thus, it is crucial to 
determine whether the VT symmetry is broken. We will see that there is a range of values of 
the parameter a in (J7J) for which the energy levels are real and this is the region of unbroken 
VT symmetry. Outside of this region some of the eigenvalues appear as complex-conjugate 
pairs. 

Let us illustrate the difference between the regions of broken and unbroken VT symmetry 
by examining some special solutions of the Schrodinger equation 

-V>"(z) + [x 6 + 2ax 4 + (4J - 1 + a 2 )x 2 }ip(x) = Eip(x), (11) 

corresponding to H in (J7J). First, consider the case J = 1. The unique eigenfunction solution 
to (fTTj) of the form in (JHJ) is ip(x) = expire 4 + ^ax 2 ) and the corresponding eigenfunction 
is E = —a. Note that E is real so long as a is real. Thus, for J = 1 there is no region of 
broken VT symmetry. 

Next, consider the case J = 2. Now, there are two eigenfunctions. The two eigenvalues 
are given by 

E = -3a± 2Va 2 - 2. (12) 

Thus, there is now an obvious transition between real eigenvalues (unbroken VT symmetry) 
and complex eigenvalues (broken VT symmetry). Evidently, the eigenvalues are real if 
a > y/2 or if a < -\/2. 

We find that for any positive integer value of J > 1, the eigenvalues E for H in (J7J) are 
entirely real if a 2 is greater than some critical value [a cr i t (J)] 2 that depends on J. These 
critical values up to J = 20 are shown in Table UJ 

Observe from TableUthat the critical values of [a cr i t ] 2 grow monotonically with increasing 
J. We have therefore also calculated the differences between successive critical values of a 2 . 
These differences also grow monotonically with increasing J, but they appear to be leveling 
off and seem to be approaching a limiting value. To see whether the differences are indeed 
approaching a limiting value as J increases, we have plotted in Fig. ^ these differences as a 
function oil/ J. This plot suggests that the differences tend to the value 12 as J — ► oo. 

To determine whether it is true that these differences really do approach limit 12, it is 
necessary to extrapolate the sequence of differences to its value at J = oo. To do so we 
have calculated the Richardson extrapolants 18] of the sequence of differences. The first 
Richardson extrapolants, Ri(J), of these differences are listed in Table ITT1 Observe that 
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TABLE I: Critical values, [a cr it(J)] 2 , of the parameter a 2 listed as a function of J. When a 2 is 
greater than this critical value, the eigenvalues of the PT-symmetric Hamiltonian H in (J7J) are 
all real. Thus, this is the region of unbroken VT symmetry. The VT symmetry is broken when 
a 2 < [a C rit(<^)] 2 - Note that the differences between successive values of [o C rit 

(J)] 2 appear to be 

approaching a limit and this is indeed the case. In fact, the numerical value of this limit is exactly 
12. Thus, for large J the critical values have the simple asymptotic behavior [a cr it (J)] 2 ~ 12J. 

the sequence Ri(J) rises more slowly and quite convincingly appears to be approaching the 
value 12. The differences Ri(J+ 1) — R\(J) between successive Richardson extrapolants are 
also shown. 

To test further the hypothesis that R\{J) tends to the limiting value 12 as J — > oo, we 
have calculated successive Richardson extrapolants of the Richardson extrapolants R\(J) in 
Table The successive extrapolants are listed in Table 11111 and they provide very strong 
numerical evidence that limj^oo ([a cr i t (J+ l)] 2 — [a cr i t (J)] 2 ) = 12. From this we conclude 
that for large J the asymptotic behavior of the critical value of a 2 is given by 

Krit(J)] 2 ~ 12J (J->oo). (13) 

Our numerical analysis provides convincing evidence that for large J the boundary be- 
tween the regions of broken and unbroken VT symmetry is given by the asymptotic behavior 
in (113)) . We will now verify this result analytically by using WKB methods [l8||. From our 
numerical analysis we know that the first eigenvalues to become complex conjugate pairs 
are always the highest, and this implies that WKB is the appropriate tool for investigating 
the VT boundary for large J. 
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FIG. 1: The differences [a cr i t (J+ l)] 2 — [a cr it(J)] 2 taken from Table [I] plotted as a function of X/J. 
Observe that as J increases, these differences tend to towards the limiting value 12. 
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TABLE II: First Richardson extrapolants Ri(J) of the sequence of differences [a cr it(J + l)] 2 — 
[a C rit(«/)] 2 taken from Tabled Notice that Ri(J) rises slowly and smoothly towards its limiting 
value 12. The differences between successive Richardson extrapolants are also listed. 

For the potential V(x) = x 6 + 2ax 4 + (a 2 + AJ — l)x 2 , the leading-order WKB quantization 
condition, valid for large n, is 

(2n + i)?r ~ J?* dx^E n - V(x) (n -> oo), (14) 

where are the turning points. Note that there is a factor of 2n + ^, rather than n + ^, on 
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TABLE III: Repeated Richardson extrapolants of the sequence of Richardson extrapolants in Table 
ILTI This table provides strong and convincing numerical evidence that Richardson extrapolants 
R\{J) tend to the limiting value 12 as J — > oo. This implies that for large J the critical values of 
a 2 grow linearly with J. See Eq. (|13|) . 



the left side of this asymptotic relation because we are counting even-parity eigenf unctions. 

For large n = J we approximate the integral in (JT%j) by making the asymptotic substi- 
tution a ~ \f~Jb, where b is a number to be determined. In order to verify the asymptotic 
behavior in (jTSj) , we must show that b = y/l2. We then make the scaling substitutions 

x = yJ 1/A and Ej ~ FJ 3/2 (15) 

because for large J we can then completely eliminate all dependence on J from the integral. 
We thus obtain the condition 

2tt= / dy^F - [y e + 2by 4 + (b 2 + 4)y% (16) 

where Ui^ — T\,iJ~ X are zeros of the algebraic equation 

y 6 + 2by 4 + (b 2 + A)y 2 - F = 0. (17) 

Next, following the analysis in Ref. Q, we assume that in this large- J limit the polynomial 
in (|T7j) factors: 

(y 2 - «) V - /3) = 0. (18) 

The correctness of this factorization assumption will be verified in the subsequent analysis. 
We then expand ()18|): 

y 6 - y\f3 + 2a) + y 2 (a 2 + 2a{3) - a 2 (3 = 0. (19) 
Comparing coefficients of like powers of y in ()17j) and (|19p. we obtain the three equations 

F = a 2 13, (20) 
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26 = -2a -P, (21) 

b 2 + 4 = a 2 + 2a(3. (22) 

Subtracting the square of Eq. (f2~Tj) from three times ((221), we § e ^ ft — at = ±V^ 2 — 12, 
and solving this equation simultaneously with (|21|) . we get expressions for a and (5: 



3a = -26 - Vb 2 - 12, (23) 



3/3 = -26 + 2V6 2 - 12. (24) 
We then substitute (J2BJ) and fl21) into (j201) to obtain 



F = -±(b - Vb 2 - 12)(26 + V6 2 - 12) 2 . (25) 

Finally, we calculate the value of the number b. Our procedure is simply to show that 
the special choice b 2 = 12 is consistent with the limiting WKB integral in Q16j) . With this 
choice we can see from f!23j) and (j2H) that a = (3 = 4/ a/3 and that (fTT)J) reduces to 

2vr = f rfy (a - y 2 ) 3 / 2 = 2 T dy(a- y 2 f' 2 . (26) 

Jy=—a J y=Q 

We simplify this integral by making the substitution y = ^/ua, and obtain 

ln = fl duu-V\l-ur/ 2 , (27) 

which is an exact identity. Thus, we may conclude that b 2 = 12. This verifies the asymptotic 
formula in (j!3)l for the location of the VT boundary. 

Furthermore, we can see that F = yV^ ~ 12.3. Thus, we obtain a formula for the 
large- J asymptotic behavior of the largest QES eigenvalue at the VT boundary: 

Ej~f^3J^ (J-oo). (28) 

The difference between this WKB calculation and that done in Ref. 0| for the Hermitian 
QES sextic Hamiltonian (JIJ is that here we have a critical value, b = y/l2, or a ~ 
This critical value defines the boundary between the regions of broken and unbroken VT 
symmetry for the PT-symmetric Hamiltonian in (jZJ. There is no analog of this boundary 
for Hermitian Hamiltonians. 
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